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Abstract Using a sample of over 70,000 red clump (RC) stars with 5-10% distance 
accuracy selected from the LAMOST Spectroscopic Survey of the Galactic Anti-center 
(LSS-GAC), we study the radial and vertical gradients of the Galactic disk(s) mainly in 
the anti-center direction, covering a significant disk volume of projected Galactocentric 
radius 7 < Rgc < 14 kpc and height from the Galactic midplane 0 < \Z\ < 3 kpc. 
Our analysis shows that both the radial and vertical metallicity gradients are negative 
across much of the disk volume probed, and exhibit significant spatial variations. Near 
the solar circle (7 < Rgc < 11.5 kpc), the radial gradient has a moderately steep, neg¬ 
ative slope of —0.08dexkpc _1 near the midplane (\Z\ < 0.1 kpc), and the slope flattens 
with increasing \Z\. In the outer disk (11.5 < Rgc < 14kpc), the radial gradients have 
an essentially constant, much less steep slope of —0.01 dexkpc 1 at all heights above 
the plane, suggesting that the outer disk may have experienced an evolution path different 
from that of the inner disk. The vertical gradients are found to flatten largely with increas¬ 
ing Rgc- However, the vertical gradient of the lower disk (0 < \Z\ < 1 kpc) is found 
to flatten with -Rgc quicker than that of the upper disk (1 < \Z\ < 3kpc). Our results 
should provide strong constraints on the theory of disk formation and evolution, as well 
as the underlying physical processes that shape the disk (e.g. gas flows, radial migration, 
internal and external perturbations). 
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1 INTRODUCTION 

The distribution of metals in the disk of a galaxy is the consequence of a variety of complicated pro¬ 
cesses, including star formation and evolution, gas flows and stellar migration, non-axisymmetric per¬ 
turbations (e.g. by the central bar, spiral arms), and interactions and accretions. The stellar metallicity 
distribution in a disk galaxy including our own Milky Way typically exhibits a (negative) gradient both 
in the radial and vertical direction. Measurements of (radial and vertical) metallicity gradients are there¬ 
fore of great importance to constrain the formation and evolution of galaxies, especially of the Milky 
Way. 
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The radial metallicity gradient of the Milky Way disk(s) has been measured using various trac¬ 
ers, including classical Cepheids (e.g. Andrievsky et al. 2002a,b,c, 2004; Luck et al. 2003, 2006, 
2011; Kovtyukh et al. 2005; Yong et al. 2006; Lemasle et al. 2007, 2008; Pedicelli et al. 2009, 2010; 
Luck&Lambert 2011; Genovali et al. 2013, 2014; Martin et al. 2015), H II regions (e.g. Deharveng et 
al. 2000; Quireza et al. 2006; Rudolph et al. 2006; Balser et al. 2011), O/B-type stars (e.g. Smartt & 
Rolleston 1997; Gummersbach et al. 1998; Rolleston 2000; Daflon & Cunha 2004; Daflon et al. 2009), 
open clusters (OCs; Friel et al. 2002,2010; Chen et al. 2003; Yong et al. 2005,2012; Carraro et al. 2007; 
Jacobson et al. 2008, 2009, 201 la,b; Magrini et al. 2009, 2010; Carrera & Pancino 2011; Frinchaboy et 
al. 2013) and planetary nebulae (PNe; Henry et al. 2004, 2010; Maciel, Lago & Costa 2005; Stanghellini 
et al. 2006; Perinotto & Morbidelli 2006; Stanghellini & Haywood 2010). The above tracers are mainly 
close to the midplane and the sample sizes are generally limited to a few hundred targets. Recently, with 
the advent of large scale spectroscopic surveys of the Milky Way, radial gradients at different heights 
(Z’s) from the midplane and vertical gradients at different projected Galactocentric radii (7?gc’s) have 
been measured with large numbers (~ tens of thousands) of field main-sequence/giant stars, such as 
main-sequence turnoff stars used by Cheng et al. (2012) from the SDSS/SEGUE (York et al. 2000; 
Yanny et al. 2009), both dwarf and giant stars used by Boeche et al. (2013,2014) from the RAVE survey 
(Steinmetz et al. 2006), red giants used by Hayden et al. (2014) from the SDSS-III/APOGEE (Allende 
Prieto et al. 2008; Ahn et al. 2014). Those tracers retain the chemical composition of the interstellar 
medium at the time of their birth. Hence, tracers of different ages probe metallicity gradients of the 
disk at different epoches. Near the solar circle (7 < Rgc < 10 kpc), the measured radial gradients of 
the disk measured with different tracers range from —0.01 to —0.09 dex kpc -1 (e.g. Cheng et al. 2012; 
Lemasle et al. 2013) and the vertical gradients range from —0.07 to —0.30 dex kpc -1 (e.g. Katz et al. 
2011; Chen et al. 2011; Kordopatis et al. 2011; Schlesinger et al. 2012, 2014). The negative radial gra¬ 
dients are generally explained with the so-called inside-out chemical evolution model (e.g, Larson et 
al. 1976; Matteucci & Francois 1989; Chiappini et al. 1997, 2001; Hou et al. 2000; Prantzos & Boissier 
2000; Molla & Diaz 2005; Fu et al. 2009) in which the inner region of the disk is chemically enriched 
faster than the outer region. 

However, there has been evidence that even close to the disk midplane, a simple linear relation may 
not be sufficient to describe the full variations of metallicity over the full range of Rgc - Using classical 
Cepheids, several studies (e.g. Andrievsky et al. 2002b; Pedicelli et al. 2009, 2010; Luck & Lambert 
2011; Genovali et al. 2013) have discovered that the radial gradient in the inner disk (4 < Rgc < 7 kpc) 
is steeper than that near the solar circle. On the other hand, a flat, nearly zero radial gradient in the inner 
disk (Rgc < 6 kpc) has been found recently by Hayden et al. (2014) using APOGEE red giants. In the 
outer disk (Rgc >10-12 kpc), most studies based on OCs find that the radial gradient tends to be close 
to zero, and is flatter than that near the solar circle. Similar flat gradients in the outskirts of the disk 
are also reported by studies using other tracers, including H II regions (Vilchez & Esteban 1996), PNe 
(Costa et al. 2004) and classical Cepheids (Andrievsky et al. 2002c, 2004; Luck et al. 2003; Lemasle 
et al. 2008). However, it is still not clear that those observed flat gradients are real or simply due to the 
coarse sampling since almost all the aforementioned studies are based on samples containing only few 
dozens of targets in the outer disk. More recently, more classical Cepheids have been found in the outer 
disk. The newly radial gradient seems to have a value comparable to that near the solar circle, lending 
further evidence that the previously reported flat gradients in the outer disk may be indeed artifact caused 
by the effects of poor sampling (e.g. Lucket al. 2011; Luck&Lambert 2011). From the theoretical point 
of view, a flat radial gradient in the outer disk can be explained by the natural outcome of disk formation 
and evolution (e.g. Carraro et al. 2007), the mixing effects due to the presence of non-axisymmetric 
structures (e.g. central bar, long-lived spiral arm; e.g. Andrievsky et al. 2004; Scarano & Lepine 2013), 
or by a merge event (e.g. Yong 2005). 

Out of the midplane, the radial gradients are found to flatten or even reverse to a positive slope as 
with the height from the midplane increases (Allende Prieto et al. 2006; Juric et al. 2008; Carrell et al. 
2012; Cheng et al. 2012; Hayden et al. 2014). The flattening of the radial gradient are possibly linked to 
the distributions of different stellar populations of the Galactic disk(s) and may possibly be explained by 
the different levels of mixture of thin- and thick-disk stars as height varies (e.g. Juric et al. 2008). The 
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vertical gradients have also been found to show spatial variations in the sense that the gradients flatten 
with increasing Rgg ( e -g- Hayden et al. 2014). 

The evolution of radial gradient is of utter importance for an understanding of the formation and 
evolution of the Galactic disk(s). However, the results of existing studies of temporal variations of the 
radial gradient of the disk are far from being unequivocal, both observationally and theoretically. Most 
studies based on OCs have found a flattening radial gradient with cosmic time by dividing the OCs into 
different age bins. The results from PNe are less clear: the flattening, steepening and unchanged are all 
reported in the literature (e.g. Maciel, Lago& Costa 2005; Stanghellini & Hay wood 2010; Henry et al. 
2010). It should be emphasized that all the tracers (OCs/PNe) used in the above analyses are younger 
than about 5-7 Gyr (see Maciel, Lago & Costa 2005). The situation of possible radial gradient evolution 
at the early epochs of disk formation is still unclear. Theoretically, different chemical evolution models 
predict different evolution of the radial gradient: flattening with cosmic time (e.g. Molla et al. 1997; 
Portinari & Chiosi 1999; Hou et al. 2000) or steepening with cosmic time (e.g. Tosi 1988; Chiappini 
et al. 1997, 2001). More recently, the radial migration has been suggested to be an important process 
in regulating the evolution of radial gradient and thus has to be included in the chemical evolution 
models (e.g. Schonrich & Binney 2009; Loebman et al. 2011; Kubryk et al. 2013; Minchev et al. 2013). 
Generally, the radial gradients of the old populations will be strongly affected and flattened by the 
migration. 

Compared to measurements of the exact slope of the metallicity gradient, studies of the spatial 
and temporal variations of the gradients are more fundamental for understanding the formation and 
evolution of the Galactic disk(s) since the variations contain information about the underlying physical 
processes that regulate the disk formation and evolution at different times, and in different regions of 
the disk. However, as described above, the currently available measurements of the spatial and temporal 
variations of the gradients are still quite controversial and in doubt. This is possibly because most of 
the previous studies have not considered the temporal variations of gradients when deriving the spatial 
variations of gradients, or, the spatial variations of gradients when deducing the temporal variations of 
gradients , given the following considerations. Firstly, the metallicity gradients in different regions of 
the disk may vary significantly, due to, for example, the effects of the non-axisymmetric perturbations 
or merge events. Secondly, as discussed earlier, different tracers probe metallicity gradients at different 
evolutionary epochs of the Galaxy. Finally, the potential biases caused by target selection effects are 
not always fully explored in many of the existing studies. This last effect is, at least partly, responsible 
for the differing results yielded by different studies. A large-scale spectroscopic survey that targets 
large numbers of stars over the full age range of the Galaxy and covers a sufficiently large volume of 
the Galactic disk(s) with sufficient sampling number density for all stellar populations concerned, is 
thus required to address the above problems properly and eventually draw a clear picture about how the 
Galactic disk(s) forms and evolves chemo-dynamically. 

The Large Sky Area Multi-Object Fiber Spectroscopic Telescopjj] (LAMOST, also named the 
Guoshoujing Telescope) Spectroscopic Survey of the Galactic Anti-center (LSS-GAC; Liu et al. 
2014; Yuan et al. 2015), as a major component of the on-going LAMOST Experiment for Galactic 
Understanding and Exploration (LEGUE; Deng et al. 2012), aims to collect low resolution (R ~ 1800), 
optical (AA3800-9000) spectra under dark and grey lunar conditions for a statistically complete 
sample of over 3 million stars of magnitudes 14.0 <r < 17.8 mag (down to 18.5 mag for limited 
fields), in a contiguous sky area of ~ 3400 sq.deg., centred on the Galactic anticentre (GAC), cover¬ 
ing Galactic longitudes 150 < l < 210° and latitudes |6| < 30°. Over 1.5 million spectra of very bright 
stars (9 < r < 14.0 mag) will also be obtained under bright lunar conditions. The survey will deliver 
spectral classifications, values of stellar radial velocity V r and stellar atmospheric parameters (effective 
temperature T e g, surface gravity log g, metallicity [Fe/H], [o/Fe ] and [C/Fe]) from the collected spectra. 
As designed, the survey will provide a large number of stars, of the order of a few millions, with a con¬ 
tiguous sky coverage of a large area of the Galactic disk(s) with a notably sampling density (hundreds 


1 LAMOST is a 4 metre quasi-meridian reflecting Schmidt telescope equipped with 4000 fibers, each of an angular diameter 
of 3.3 arcsec projected on the sky, distributed in a circular field of view of 5° in diameter (Cui et al. 2012). 
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of stars per squared degree) yet using a simple but nonetheless nontrivial target selection algorithms 
(random selection in the colour-magnitude diagram, i.e. r versus g — r and r versus r — i). The survey 
will thus enable simultaneous measurements of the radial and vertical metallicity gradients across a 
large volume of the Galactic disk(s). The LSS-GAC regular survey was initiated in September 2012. It 
is expected to last for five years. 

Based on the large stellar spectroscopic samples obtained by the LSS-GAC, we have started a series 
of investigations to systematically study the spatial and temporal variations of metallicity gradients of the 
Galactic disk(s). Xiang et al. (this volume, hereafter Paper I) examine the spatial and temporal variations 
(in both radial and vertical directions) of metallicity gradients for 7.5 < Rgc < 13.5 kpc and \Z\ < 
2.5 kpc using a sample of nearly 0.3 million main-sequence turnoff (MSTO) stars of ages ranging from 2 
to over 11 Gyr, selected from the LSS-GAC DR2 (Xiang et al. 2015a, in preparation). The current study, 
the second one of this series, selects over 70, 000 red clump (RC) giants of intermediate- to old-ages 
also from the LSS-GAC DR2. The bright intrinsic luminosities of RC stars allow us to investigate the 
(radial and vertical) gradients over a slightly larger volume of the disk, covering 7.0 < Rgc < 14.0 kpc 
and \Z\ <3.0kpc, than possible with the MSTO star sample used by PaperI. The current sample is 
thus particularly suitable for the study of the outer disk and provides a unique opportunity to address 
the question whether the radial gradient of the outer disk is flatter than that near the solar circle as 
described above. The paper is organized as follows. We describe in detail the selection of RC sample 
in Section 2. We present our results in Section 3, followed by discussions in Section4. Conclusions are 
given in Section 5. 

2 DATA AND SAMPLE SELECTION 

2.1 Data 

The data used in the current study are mainly collected during the Pilot (Sep. 2011-Jun. 2012 ) and 
first-two-year (Sept. 2012-Jun. 2014) Regular Surveys of LSS-GAC. The stellar radial velocity V r and 
stellar atmospheric parameters (T e ff, log g, [Fe/H]) are determined with the LAMOST Stellar Parameter 
Pipeline at Peking University (LSP3; Xiang et al. 2015b,c), and have an overall accuracy of 5-10kms _1 , 
150K, 0.25 dex, 0.15 dex, respectively, for a spectral signal-to-noise ratio per pixel (~ 1.07 A) at 4650 A, 
S/N (4650 A) >10. Stellar parameters, together with values of interstellar extinction and stellar distance 
estimated with a variety of techniques of over 0.7 million stars observed in the Pilot and first-year 
Regular Surveys are now publicly available as the first data release of value-added catalog of LSS-GAC 
(LSS-GAC DR1; Yuan et al. 2015). The data of another 0.7 million stars with S/N (4650 A) > 10 from 
the second-year Regular Survey will soon be publicly available as the second release of value-added 
catalog of LSS-GAC (Xiang et al. 2015a, in preparation). In addition to the LSS-GAC data, the current 
study has also included about 80,000 stars from the spheroid parts (|6| > 20°) of the LEGUE and 
50,000 stars from LAMOST- Kepler Field Survey (De Cat et al. 2014). The latter are particularly 
useful for the RC selections (see Section 2.2). The additional spectra from the latter two surveys are also 
processed in the same way as for the LSS-GAC sources with the LSP3. 

2.2 The RC sample 

2.2.1 An overview of RC stars 

The helium burning RC star^ are an easily recognizable feature in the colour-magnitude diagram 
(CMD) of intermediate- to old-age, low-mass but metal-rich stellar populations. RC stars are consid¬ 
ered as standard candles since their absolute luminosities are fairly independent of the stellar chemical 
composition and age (e.g. Cannon 1970; Paczynski & Stanek 1998). With the accurate calibration of 
absolute magnitudes of several hundred local RC stars using the Hipparcos (ESA 1997) parallaxes be- 

- RC stars of interest here refer to those in which the helium burning is ignited following a helium flash that breaks the 
degeneracy of the helium core. 
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coming available, RC stars have become popular distance indicators and been widely used to obtain 
precise distances to the Galactic Center (Paczynski & Stanek 1998) and the Local Group of galaxies 
(e.g. LMC by Laney et al. 2012; M31 by Stanek & Granvich 1998). The key of this method is a notable 
over-density feature easily identifiable in the CMD of stars in those distant objects due to the presence 
of a large number of RC stars at nearly the same distance. 

However, identifying individual RC stars amongst the numerous field stars is not an easy task. Field 
RC stars spread over a wide range of distances and there is no apparent over-density in the CMD of 
field stars that can be used to single out the RC members. Fortunately, high-resolution spectroscopic 
analysis of nearby RC stars show that they distribute in a relative “small box” in the T e g - log g diagram 
(i.e. HR diagram) of 4800 < T e g < 5200 K and 2.0 < log g < 3.0 (e.g. Puzeras et al. 2010). This is sim¬ 
ply because T e g is an excellent proxy of stellar color whereas log g is sensitive to the stellar absolute 
luminosity. Recently, with values of T e g and log g available from large-scale spectroscopic surveys for 
large numbers of field stars, it has become feasible to select a large number of field RC candidates. RC 
candidates thus selected have been widely used to study the metallicity gradients and stellar kinematics 
of the Galactic disk(s) (e.g. Siebert et al. 2011; Bilir et al. 2012; Williams et al. 2013; Bienayme et 
al. 2014). However, stars falling inside that “small box” are not purely RC stars and have significant 
(~ 60%; Williams et al. 2013) contamination from red giant branch (RGB). The differences in absolute 
magnitudes between RC and RGB stars can be larger than 1 mag and this can lead to large systematic 
errors in distances of the selected RC sample stars. More recently, a new method has been proposed 
by Bovy et al. (2014; hereafter B14) to select a clean RC sample from large-scale spectroscopic data. 
The method first separates the RC-like and RGB stars using cuts on a metallicity-dependent T e g - log g 
diagram, assisted by theoretical using stellar isochrones and calibrated using Kepler high quality aster- 
oseismic log g (e.g. Creevey et al. 2013) as well as high precision stellar atmospheric parameters from 
the APOGEE survey (Pinsonneault et al. 2014). Secondly, secondary RC stars (i.e. high mass helium 
burning stars ignited nondegeneratedly) are then removed from the RC-like stars via cuts on a metallic¬ 
ity (Z) -color (J — K s )o diagram. The expected purity of the final RC sample is > 93 % and the typical 
distance errors are within 5 - 10 %. The key point of this new method in selecting a clean RC sample is 
accurate values of surface gravity logg. Unfortunately, it is difficult to apply this method to the current 
LSS-GAC data given the relatively large systematic plus random errors (~ 0.2-0.4dex) of log g esti¬ 
mates delivered by the LSP3, as shown by a recent comparison of LSP3 and asteroseismic log g values 
for common objects in the LAMOST -Kepler fields (Ren et al. 2015, in preparation). As we shall show 
later (Section2.2.2; Fig. 4) there are also significant deviations between T e g-log g values yielded by the 
LSP3 and those predicted by the PARSEC stellar isochrones (Bressan et al. 2012) of LSS-GAC DR2 
red giant stars for a typical metallicity [Fe/H] = —0.3. Therefore, for selecting a clean RC sample from 
the LSS-GAC data set, it is necessary to improve the accuracy of LSP3 log g estimates.. 

2.2.2 Re-determinations ofLSP3 log g estimates with the KPCA algorithm 

Asteroseismology is very powerful for determination of stellar surface gravity log g, and is capable of 
yielding values much more accurate than possible with the traditional spectroscopic method. Thanks to 
the Kepler mission (Borucki et al. 2010), asteroseismic values of log g of tens of thousands stars are 
now available, with a typical error of ~0.02dex (Huber et al. 2014). This number of stars is however 
still orders of magnitude smaller than those already been targeted by large-scale spectroscopic surveys, 
including LSS-GAC. It will be thus very useful to calibrate and improve the accuracy of spectroscopic 
estimates of log g using the much more accurate asteroseismic values. Very recently, Liu et al. (2014) 
improve the accuracy (by a factor of two) of log g estimates of LAMOST Stellar Parameter Pipeline 
(LASP; Wu et al. 2011; Luo et al. 2015) using a support vector regression model trained by astero¬ 
seismic log g measurements of thousands of giants from the LAMOST-Aep(er fields. Here, a similar 
effort is presented for the LSP3 log g estimates . We re-determine and improve the accuracy of LSP3 
estimates of logg based on a Kernel Principal Component Analysis (KPCA; Scholkopf et al. 1998) 
method trained with thousands of giants in the LAMOST-A epler fields with accurate asteroseismic 
log g measurements. 
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Fig-1 Distributions of values of residual, log <7kpca— log <7ast- of the training sample as 
a function of LSP3 T e s and [Fe/H] (Panels a and b) and of asteroseismic log g (Panel c). 
Panel (d) shows a histogram distribution of the residuals (black line). Also overplotted in red 
is a Gaussian fit to the distribution. The mean g = 0.00 dex and dispersion cr = O.lOdex of 
the fit, as well as the number of stars used, N = 1355, are marked. For the current KPCA 
analysis, Npc = 25 is assumed. See text for detail. 


The Principal Component Analysis (PCA) is a classical statistical method that can convert the ob¬ 
servations (e.g. spectra) into a set of linearly uncorrelated orthogonal variables or principal components 
(PCs). The method has widely been used to estimate stellar atmospheric parameters (e.g. Singh et al. 
1998; Rees et al. 2000; Re Fiorentin et al. 2007). The KPCA remains the concept of PCA but can be 
extended to nonlinear feature extraction using the kernel techniques. Considering the nonlinear depen¬ 
dence of spectral features on the log g, the KPCA is more suitable than the standard PCA method for 
the case of log g estimation. 

To apply the KPCA method to the log g determinations for a whole spectroscopic sample, it is 
necessary to construct a tight relation between the LAMOST spectra (at the present, only the blue- 
arm spectnO between 3930 < A < 5500 Aare used) and the corresponding asteroseismic values of log g 
with training data. To build up a training sample, we have cross matched our sub-sample of about 
the 50,000 stars in LAMOST -Kepler fields with the currently available largest asteroseismic log g 
(hereafter log <?ast) sample (~ 16,000 stars) from Huber et al. (2014). In total 3562 common sources 
of spectral S/N (4650 A) > 10 are identified. For self-consistency, we have used the LSP3 estimates 
of T e ff and the observed frequencies of maximum power (I'm ax) from Huber et al. (2014) to re-estimate 
values of log g for the 3562 common sources, using the asteroseismic scaling relation !/ max oc g T e g 0 ' 3 * 5 
(Brown et al. 1991; Belkacem et al. 2011). We then trim this sample of the common sources requiring 
that LSP3 log g < 3.8, 3500 < T e s < 6000K and [Fe/H] > —1.0, for the following three reasons: 1) The 
current study concerns mainly the estimation of log g for red giant stars including RGB and RC stars; 
2) Currently only a small number of dwarfs have asteroseismic log g determinations; and finally, 3) The 
surface gravity of metal-poor stars with [Fe/H] < —1.0 are poorly determined currently by astroseismic 
scaling relation (Epstein et al. 2014). The trimming leaves 3267 sources. The 3267 sources are then 
further divided into two groups: a training sample and a control sample. For the training sample, we 
select 1355 stars with spectral S/N (4650 A) > 50 to minimize the effects of spectral noises and of errors 
of LSP3 stellar atmospheric parameters T e g and [Fe/H] in the training process. The remaining 1912 
stars constitute the control sample, which will be used to test the accuracy of log g determined with 
KPCA algorithm. 

To extract information of the stellar surface gravity from the LAMOST spectra, we assume a KPCA 
model (with a Gaussian kernel) of 25 Principle Components (PCs), iVpc = 25, when constructing the 

3 In principle, red-arm spectra should also be included given that the stars of concern here are mainly red giants that are 
bright in the red and thus generally have high spectral SNRs. However, considering the sky subtraction of the current LAMOST 

2D pipeline (Luo et al. 2015) which still suffers from some problems in the red and the fact that most stars typically have few 

spectral features in the red, we have opted to consider the blue-arm spectra in the current study. In the future, with improved sky 
subtraction, we expect including the red-arm spectra can help improve the results. 
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Fig. 2 Differences of KPCA log g estimates and asteroseismic values for the control sample, 
plotted as a function of spectral S/N (4650 A). The blue dashed lines indicate differences 
A log g = ± 0.1 dex. The means and standard deviations of the differences in the individual 
S/N (4650 A) bins (binsize = 8 ) are oveplotted with dots and error bars, respectively. Again, 
Npc = 25 is assumed. 
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Fig. 3 Differences of KPCA log g estimates yielded by duplicate observations of similar 
S/N (4650 A) plotted as a function of S/N (4650 A). The blue dash lines indicate differences 
A log g = ± 0.1 dex. The means and standard deviations of the differences in the individual 
S/N (4650 A) bin (binsize = 8) are oveplotted as red dots and error bars, respectively. Again, 
Npc = 25. 
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Fig. 4 Pseudo-color HR diagrams of red giant stars from the LSS-GAC DR2. Only stars of 
LSP3 metallicity [Fe/H] between —0.4 and —0.2 (where the metallicity distribution of our red 
giant sample peaks) are included. Effective temperatures are those given by the LSP3. In the 
left panel, log g values are from the LSP3 , whereas those in the right panel are newly derived 
with the KPCA method. Also overplotted are the PARSEC stellar isochrones for [Fe/H] = 
—0.3. Different colors represent different stellar ages (1-2 Gyr, green crosses; 2-6 Gyr, blue 
crosses; 6-10 Gyr, red crosses). 


relation between log c/ast and the LAMOST blue-arm spectrum with the training sample. Ape = 25 is 
chosen in the KPCA model for a tradeoff between reducing the training residual and retaining a relative 
high accuracy of log g estimate for LAMOST spectra at low spectral S/N’s (see Appendix A). In the 
training process, T e g and [Fe/H] are fixed to the value yielded by the LSP3. The model fits the data well, 
yielding residuals of only ~ 0.1 dex, as shown in Fig. 1. In addition, the Figure shows that the residuals 
exhibit no obvious systematic trend as the LSP3 T e g and [Fe/H] vary. But we caution that the residuals 
show a weak trend as the asteroseismic log g varies from 2.5 to 3.5. The reason is unclear but the trend 
do not affect the selection of RC stars (see Fig. 5). To test the accuracy of KPCA log g values, we apply 
the relation yielded by the training sample to the control sample. The results are shown in Fig. 2. The 
differences between KPCA and asteroseismic log g are less than 0.15 dex for spectral S/N (4650 A) > 18 
and less than 0.10 dex for S/N (4650 A) > 30. Thus the accuracy of KPCA log g estimates are about a 
factor of two better than the LSP3 values. The test shows that there is no obvious systematics over the 
whole parameter space (atmospheric parameters as well as S/N ratios) spanned by the control sample. 

We now apply the above model to the LAMOST blue-arm spectra and re-determine log <7 val¬ 
ues for all LSS-GAC DR2 red giant^ using the relation constructed from the training sample. This 
yields newly derived KPCA log <7 values for over 0.35 million LSS-GAC DR2 red giantfQ with spec- 


4 Red giants are selected from the whole LSS-GAC DR2 sample using the same set of stellar atmospheric parameter cuts as 
for the training sample, i.e. log g < 3.8, 3500 < T e ff < 6000 K and [Fe/H] > —1.5 . 

5 This includes targets from the LSS-GAC as well as stars from the spheroid part of the LEGUE and the LAMOST -Kepler 
fields. 
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Fig. 5 Distribution of stars on the T e g - log g HR diagram predicted by the PARSEC stellar 
evolution model with assumptions as described in the text (grey scale), compared to that of 
the LE sample stars (colored points), for two metallicity bins as marked. Values of T e s and 
[Fe/H] of the LE sample stars are from LSP3 and those of log g from the KPCA method. 
Red, magenta and blue points represent stars of RC, RGB and unknown evolutionary stages, 
as classified by Stello et al. (2013), respectively. The blue dashed lines represent the cuts that 
separate RC stars and the less luminous RGB stars as given in Eqs. (1) and (2). The green 
dashed lines give the cuts adopted by B14. 


tral S/N (4650 A) > 10. As another test of the accuracy of newly derived KPCA log g estimates, we 
plot the differences of KPCA log g values yielded by ~ 6500 pairs of duplicate observations of similar 
S/N (4650 A) as a function of S/N. The results are quite similar to those found for the control sample. 
Again, this indicates that our newly derived KPCA log g values are about a factor of 2 more accurate 
than those of LSP3. We also find that the HR diagram constructed using the LSP3 T e g and the newly 
derived KPCA log g for our sample stars is in excellent agreement with the prediction of the PARSEC 
stellar evolution model, as shown in Fig. 4. As mentioned before, the position of RC clump in the origi¬ 
nal HR diagram constructed with LSP3 T e g and log g deviates significantly from that predicted by the 
PARSEC model. 

2.2.3 Selection ofRC stars 

As described above, the newly derived KPCA log g values, with much improved accuracy (both in ran¬ 
dom and systematic), have much reduced the differences between the observed and predicted position 
and morphology of the RC clump on the HR diagram, thus enabling us to select a clean RC sample 
following the method proposed by B14. Firstly, to obtain a set of cuts that distinguish RC and RGB 
stars on the T e ff-logg diagram, we cross-identify red giants from the LAMOST -Kepler fields with 
those from Stello et al. (2013), in which the evolutionary stages of stars have been classified based on 
the frequencies and spacings of gravity-mode measured from the Kepler light curves. This leads to 
common 600 stars with LAMOST spectral S/N(4650A) > 40 (hereafter the LAMOST-EVOLUTION 
sample, the LE sample for short). The S/N cut is to ensure that the atmospheric stellar parameters are 
accurate enough, especially the KPCA log g values (better than 0.1 dex under this cut). The LE sample 
stars are mostly classified by Stello et al. (2013) as either ‘RC’ or ‘RGB’, only a few as ‘UNKNOWN’. 
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(J-K .) 0 


Fig. 6 Distribution of RC-like stars in the color (J — K' s )o Z plane, using 2MASS photo¬ 
metric measurements and LSP3 metallicities. Dots in red, blue and green represent main RC, 
secondary RC and RGB stars, respectively. The red dashed lines are developed by B14 to 
eliminate the contamination of secondary RC stars. The blue dashed line is our revised cut 
developed in the current study in replacement of the lower red dashed line. 


Unlike B14, who use cuts defined with the asteroseismic log g values, we opt to define the cuts using 
the LE sample, relying directly on log g values determined with the KPCA method, together with T e g 
and [Fe/H] yielded by the LSP3. With this change, there is no need to introduce an additional cut that 
correct for the systematic difference between the asteroseismic log g [e.g. the additional cut given by 
Eq. (9) of B14], nor the need to account for the possible systematic differences between the PARSEC 
and spectroscopic values of the other atmospheric parameters ( T e g and [Fe/H]). The cuts defined here 
are therefore more straightforward and appropriate for real observed data. Fig. 5 compares the distribu¬ 
tion on the 7’ e ff - log g HR diagram of the LE sample stars and that predicted by the PARSEC stellar 
isochrones for two metallicity bins. The PARSEC distribution is generated assuming, similar to B14, a 
lognormal Chabrier (2001) for initial mass function (IMF), a constant SFH with ages less than 10 Gyr 
and metallicity distribution functions (MDFs) that match with that of the LE sample stars in the two 
metallicity bins. As Fig. 5 shows, the distribution of RC-like stars, with log g values roughly between 
2.5 and 2.9, and that of RGB stars, as classified by Stello et al. (2013), is roughly in agreement with 
the predicted distribution from the PARSEC stellar evolution model except for some small systematic 
differences between the LSP3 and PARSEC values of T e g, and is well separated from that of RGBs 
as well. Finally, ignoring those stars classified as ‘UNKNOWN’ evolutionary stage, we derive the fol¬ 
lowing cuts as given by the following two equations that maximize the product of completeness f comp 
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Table 1 Values of purity and completeness of selection of RC-like stars under different un¬ 
certainties errors of log g measurements for the cuts given by Eqs. (1) and (2). 


trios q ( dex ) 

on 

013 

030 

035 

030 

S/N“ 

> 40 

20 

18 

15 

10 


94 ± 1% 

89 ± 1% 

83 ± 2% 

78 ± 2% 

74 ± 2% 

/comp 

96 ± 1% 

90 ± 2% 

83 ± 2% 

78 ± 2% 

74 ± 3% 


a S/N ratios required by the KPCA method to achieve the corresponding 
uncertainties of log g measurements listed in the first row of the Table. 


(= N^/N^ ~ 96 %) and purity / purity (= 1 - ~ 95 %), 

1-8 < logg < 0.0009dexK" 1 {T eff -T e ] | ef ([Fe/H])} + 2.5, (1) 

where 

r eff ef ([Fe/H]) = -876.8 K dex " 1 [Fe/H] + 4431IC . (2) 

represents the total number of real RC stars, and and /V B q B the numbers of RC and RGB stars 

that pass the cuts of Eqs. (1) and (2), respectively. J is the sum of and N B q B , that represents 
the total number of stars passing the cuts. 

Fig. 5 shows that the cuts developed by B14 (green lines in the Figure) is not fully suitable for our 
data because their cuts reject some real RC stars in the lower temperature range but include some RGB 
stars as RC stars in the higher temperature range. We have also plotted Fig. 5 using the asteroseismic 
values of log g from Huber et al. (2014) for the FE sample stars do not find any significant difference 
from what shown here. This indicates the difference between our current cuts and those of B14 can not 
be caused by the usage of KPCA log g. This difference in cuts is possibly caused by the systematic 
differences in values of T e g as yielded by the FSP3 and given by PARSEC, as mentioned earlier. 

To examine the performance of our cuts for different errors of log g measurements, we use a Monte 
Carlo method to simulate the selection process under different errors of log g and derive the correspond¬ 
ing values of / comp and / pU rity- Here we assume that the uncertainties of the FE sample stars have a 
minimum value of 0.1 dex (almost the best accuracy achievable using the KPCA method), based on the 
test results presented in Section 2.2.2 using the control sample and duplicate observations. Values of 
/comp and /p Ur ity for larger errors are then derived by assigning larger errors randomly to log g values 
of the FE stars and repeating the selection with Eqs. (1) and (2) by 10,000 times. The results are pre¬ 
sented in Table 1. In Table 1, we have also listed the spectral S/N(4650 A) ratios required by the KPCA 
method in order to achieve the listed uncertainties of log g measurements. With the cuts developed in 
the current work, the amount of contamination from RGB stars is twice less than that for selection using 
a “small box” such as that adopted by Williams et al. (2013), even for log g measurement uncertainties 
as large as 0.3 dex , an accuracy achievable with the KPCA method even for a spectral S/N (4650 A) 
ratio as low as 10 . 

As shown in Fig. 5, the region where RC-like stars fall show two distinct features: a clump of stars 
of almost constant log g = 2.5 and 4600 < T e g < 4850 K (i.e. the main RC stars) and a skew tail of stars 
with 2.5 < log g < 2.9 and 4800 < T e g < 5000 K (i.e. the secondary RC stars of fainter luminosities). 
The secondary RC stars are not standard candles and their luminosities are a function of mass (age). We 
need to remove them from our RC-like star sample to ensure that the distances of the sample stars can 
be determined precisely and simply. 

To remove secondary RC stars and low surface gravity RGB stars, we adopt the color ( J — K s )q - 
metallicity 20 cuts as developed by B14 using the PARSEC stellar evolution model. The cuts are aimed 
to select stars that have a nearly constant absolute magnitude with 1 a dispersion < 0.1 mag and elim¬ 
inate those non-standard candles, including secondary RC, RGB and AGB stars. For a more detailed 
description, please refer to the Section 2.2 of B14. Since the cuts are purely based on stellar evolution 
model, empirical examination is required to eliminate the potential discrepancies, if any, between the 

6 Z is converted from [Fe/H] using the equation given by Bertelli et al. (1994), assuming Z q = 0.017. 
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model and observations. To test the cuts, we again use the LE sample. For consistency, we first apply 
the cuts of B14 on the T e g - log g diagram to the LE sample to select RC-like stars using the asteroseis- 
mic log g values. After corrected for the interstellar extinction as estimated with the ‘star pair’ method 
(see Yuan et al. 2015), the distribution of selected RC-like stars in the (J — K s )q-Z plane is shown 
in Fig. 6 . The photometric magnitudes of the stars are taken from 2MASS (Skrutskie et al. 2006) and 
their metallicities from the LSP3. By definition, secondary RC stars have masses > 2M e and main 
RC stars have masses < 2 Mq. Again with masses estimated also by Stello et al. (2013), we divide the 
RC-like stars into two groups: main and secondary RC stars. As Fig. 6 shows, the contaminations of 
RGB stars, as expected, is negligible but that of secondary RC stars is significant. Luckily, compared to 
main RC stars, secondary RC stars belong to a young and metal-rich population of large masses As a 
consequence, they occupy the bluer and more metal-rich parts of color-metallicity diagram. Main and 
secondary RC stars can be well separated using a cut given by the following equation as proposed by 
B14 (the upper red dashed line in Fig. 6 ), 

X < 2.58[(J - K s )o - 0.400] 3 + 0.0034. (3) 

However, Fig. 6 also shows that the lower cut of B14 (the lower red dash line in Fig. 6 ), designed to 
reject high surface gravity RGB stars, seems to also eliminate a considerable fraction of real main RC 
stars. This is possibly caused by discrepancies between model and observations. To correct for this, we 
adjust the lower cut of B14 to (the blue dashed line in Fig. 6 ), 

X > 1.21[( J - K s ) 0 - 0.085] 7 8 9 + 0.0011. (4) 


2.2.4 The RC sample 

With 7’eff, [Fe/H] estimated with the LSP3, log g derived with the KPCA method and (J — K s )o cal¬ 
culated from 2MASS photomtrjQ after corrected for extinction as estimated with the ‘star pair’ method 
(Yuan et al. 2015), we apply the cuts of Eqs. (l)-(4) to the LSS-GAC DR2 red giant sample and ob¬ 
tain a clean RC sample of over 0.11 million stars of S/N (4650 A) > 10. The distances of these RC 
stars are derived using the recent calibration, Mk s = —1.61 mag, for nearby RC sample (Laney et 
al. 20120 Since the intrinsic scatter of absolute magnitudes of RC stars is estimated to be within 
0.1 mag, the distances derived are expected to have uncertainties no more than 5-10%, given a typ¬ 
ical photometric error of ^0.05 mag in iC s -band and an extinction error of ~ 0.04 mag in E(B V) 
(see Yuan et al. 2015). Proper motions of the sample stars are taken from the UCAC4 (Zacharias et 
al. 2013) and PPMXL (Roester et al. 2010) catalogue.s Given the high distance precision and large 
spatial coverage (Fig. 7) of the current RC star sample, it is not only useful to address the problem 
of metallicity gradients of the Galactic disk of concern here, but can also tackle a variety of problems 
with regard to the Galactic structures and dynamics. The sample will be publicly available at website: 
(http://162.105.156.2 4 9/site/RC_Sample). 

In order to obtain accurate metallicity gradients of the Galactic disk(s) of concern here, we restrict 
the current analysis to a subset of the RC sample stars with S/N (4650 A) > 15, to ensure a sufficiently 
high accuracy of [Fe/H] determinations, and to eliminate possible contaminations from halo stars of 
with [Fe/H] < —1.0 dex (also the lowest applicant range in [Fe/H] of KPCA method in log g estimates). 
Finally, a total of 74,280 unique RC stars are selected. The number density distributions of this adopted 
sample of RC stars in the X-Y and X-Z Galactic planes are presented in Fig. 7. Here X, Y and 
Z are coordinates of a right-handed Cartesian coordinate system with an origin at the Galactic center. 
The X-axis passes through the Sun and points towards the Galactic Center, the Y-axis is in the di¬ 
rection of Galactic rotation and the X-axis points towards the North Galactic Pole. The Sun is located 
at (X , Y, Z) = (—Rq, 0, 0), where Rq = 8 kpc. As Fig. 7 shows, our RC sample covers a volume of 
—16 < X < —6 kpc, —4 < Y < 5 kpc and —3 < Z < 3 kpc. 

7 Only stars with pli_qual flagged as ‘A’ in both J and K s bands are included. 

8 B14 assign absolute magnitudes to RC stars based on their positions in the (J — K s )o — Z plane. However, we find that the 

absolute magnitudes of RC stars predicted by the stellar evolution model are almost constant except for those super metal-rich 
ones (Z > 0.04, i.e. [Fe/H] > 0.4 - 0.5 dex; see Fig. 3 of B14), which are very rare in our RC sample. 
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Fig. 7 Greys cale number density distribution of RC sample stars in the X — Y (left panel) and 
X-Z (right panel) planes. The Sun is located at (X, Y, Z ) = (—8.0, 0.0, 0.0) kpc. The stars 
are binned by 0.25x0.25 kpc 2 in both diagrams. The densities are shown on a logarithmic 
scale. Histogram distributions along the X, Y and Z axes are also plotted. 


We note that the current RC sample is largely a magnitude limited sample. Considering that the 
absolute magnitudes of RC stars are quite insensitive to both metallicity and age, we believe that any 
potential population effects that may be present in our sample are unlikely to significantly affect the 
metallicity gradients derived from this sample. In addition, the cuts to define the RC sample in the 
current selection excludes stars of ages younger than 1 Gyr (i.e. those relatively massive and young 
secondary RC stars) and halo stars of metallicities poor than — l.Odex, the current sample is expected 
to be mainly composed of stars of intermediate-age stars with typical ages between 1.0 and 4.0 Gyr if 
one assumes a flat SFH and a solar-neighborhood MDF (Casagrande et al. 2011; see Section 5 of B14 
for a more detailed discussion). 

3 RESULTS 

3.1 Spatial distribution of mean metallicities 

The spatial distribution of mean metallicities of the RC sample in 0.40 x 0.25 kpc -2 bins, with at least 
five stars in each bin, in the projected Galactocentric radius ( Rgc ) and the distance from the midplane 
i\Z\) is presented in Fig. 8. For each bin, the mean [Fe/H] is derived. Fig. 8 shows several obvious 
features: 

1) A gradient in both the radial and vertical direction is evident, especially near the solar circle; 

2) By the midplane, the radial gradient near the solar circle is steeper than that in the outer disk (see 
Section 3.2); 

3) The radial gradient flattens as one moves away from the midplane (see Section 3.2); 

4) The vertical gradient near the solar circle is steeper than that at large lice (see Section 3.3). 

To describe those features of the spatial metallicity distribution in a quantitive manner, we mea¬ 
sure the radial and vertical gradients separately for the different regions of the disk in the following 
subsections. 

3.2 Radial metallicity gradients 

The radial and vertical gradients, A[Fe/H]/ ARqc and A[Fe/H]/AZ, are generally correlated with 
each other. Fortunately, with a large, contiguous disk volume sampled by a large number of stars in 
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Fig. 8 Mean metallicity [Fe/H] as a function of Rgc and \Z\ deduced from our RC sample, 
binned by 0.40 x 0.25 kpc 2 in Rgc and \Z\, respectively. 


hand as in the current case, it is now possible to determine the radial (vertical) gradient in thin vertical 
(radial) slice (annulus), and thus quantify the variations in both directions independently. 

To derive the radial gradient from our RC sample, we divide the stars into bins of constant height, 
A\Z\ = 0.2 kpc, except for the inner most bin of \Z\ <0.1 kpc and the two outer most bins of 1.9 < 
\Z\ < 2.4kpc and 2.4 < \Z\ < 3.0, respectively. For each bin, a straight line is adopted to fit the 
mean metallicities of the individual radial annuli as a function of Rgc- The slope of the straight line 
gives the radial gradient. To ensure a sufficient number of stars in all radial annuli, the binsize in radial 
direction is allowed to vary. We require that the radial binsizes are no thinner than 0.2 kpc and each 
radial annulus contains no less than 50 stars. 

Motivated by the previous findings that the metallicity variations in the radial direction across the 
entire range of the disk may not be captured by a single slope, we, inspired by Magrini et al. (2009), 
divide the whole disk into three parts: Region A of 3 < Rgc < 6 kpc, RegionB of 6 < Rgc < 
11.5 kpc and Region C of 11.5 < Rgc < 22 kpc. Region A is not covered by our RC sample, thus 
we only determine the radial gradients of Regions B and C for the aforedefined individual height bins. 
The results are presented in Fig. 9. The Figure shows that, close to the midplane {\Z\ < 1.0 kpc), 
the gradients of Region C are flatter than those of Region B. For the upper disk (\Z\ > 1.0 kpc), 
the radial gradients found for Regions B and C become comparable. Fig. 10 plots the radial gradients 
measured Regions B and C as a function of \Z\. For Region B, the radial gradient flattens with increasing 
\Z\. It has a negative gradient of —0.09 kpc dex -1 at midplane and a marginally positive of gradient 
0.02 kpc dex -1 at ~ 2.5 kpc away from the midplane. For Region C, the gradient is essentially constant, 
value of—0.014 ± 0.008dexkpc _1 . 

3.3 Vertical metallicity gradients 

To study the spatial variations of vertical gradients and also to minimize the effects of radial and vertical 
gradients on the determinations of vertical gradient, we divide the RC sample into annuli of width 1 
kpc in the Rgc direction, except for the most distant annulus of 13 < Rgc < 15 kpc. To derive the 
gradients, we calculate the mean metallicities of the individual height bins. The binsize in height \Z\ is 
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Fig. 9 Metallicity distribution as a function of Rgc in the individual disk height bins derived 
from our RC sample. The red dashed lines divide the disk into two parts: Region B near the 
solar circle and region C for the outer disk, as marked on top of the two columns of plots. Red 
boxes are the mean metallicities in the individual radial annuli. The metallicity gradients of 
Regions B and C are fitted separately, each by a straight line. The results are shown by red 
lines. The gradient s and errors thus derived are also marked near the bottom of each panel. 


allowed to vary but set to be no less than 0.2 kpc and each bin contain s no less than 50 stars. Again, a 
simple straight line is used to fit the data, and the gradient is given by the slope of the fit. 

The results are presented in Fig. 11. The derived vertical gradients as a function of Rgc are pre¬ 
sented in Fig. 12. It shows that the vertical gradient flattens with increasing Rgc- To better understand 
the spatial variations of vertical gradient, especially if there is any difference between the thin and 
thick disks, we further calculate the vertical gradients for the lower ( \Z\ < lkpc) and upper disk 
(1 < \Z\ < 3kpc) separately. The separation is chosen to be 1 kpc, where is the transition zone 
between the thin and thick disks given their scale heights (300pc for the thin disk and 900 pc for the 
thick disk; Juric et al. 2008). As Fig. 15 shows, we find that the vertical gradient of the lower disk flat¬ 
tens quicker than that of the upper disk. At Rgc > 10.5 kpc, the vertical gradients of the lower disk 
diminish to zero and even become positive, while those of the upper disk remain negative and change 
only slowly with Rgc- 

4 DISCUSSIONS 

The LSS-GAC DR2 RC sample presented here allows us to investigate the radial and vertical gradients 
over a large volume of the Galactic disk, especially of its outer region. In the following subsections, we 
will compare our results to previous determinations and discuss their implications for the Galactic disk 
formation and evolution. 
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Fig. 10 Radial gradients, A[Fe/H]/Al?GC> as a function of \Z\ yielded by the RC sample 
for disk Regions B (red dots) and C (blue dots). The lines connecting the dots have been 
smoothed over 3 adjacent points. The symbols in grey represent measurements from the lit¬ 
erature, all derived from tracers near the solar circle (i.e. Region B of this study). The star is 
the measurement by B14 using APOGEE RC stars, triangles the measurements by Anders et 
al.(2014) using APOGEE giants and boxes the measurements by Boeche et al. (2014) using 
RAVE RC stars. 


4.1 Radial metallicity gradients 

The gradients near the solar circle (i.e. Region B) are well determined in the current study given the 
large number of RC stars with accurate distances assembled here. Near the midplane (\Z\ < 0.1 kpc), 
the radial gradient measured from our RC sample has a value of 0.082 ± 0.003 dex kpc -1 , which is in 
excellent agreement with the result of 0.09 ± 0.01 dex kpc -1 obtained by B14, also using RC stars near 
the midplane ( \Z\ < 0.05 kpc) with a similar coverage of Rqc- Our results show that the radial gradient 
flattens with increasing disk height \Z\. A similar trend has been found by a number of recent studies 
(e.g. Cheng et al. 2012; Anders et al. 2014; Boeche et al. 2014; Hayden et al. 2014). In Fig. 10 we have 
also overplotted the measured radial gradients measured at different heights \Z\ by Anders et al. (2014) 
using APOGEE red giants and by Boeche et al. (2014) using RAVE RC stars. Their measurements are 
generally in good agreement with our results. If one assumes that thin disk stars concentrate more closely 
to the midplane than thick disk stars, the trend of variations of radial gradient with \Z\ may suggest that 
the thin disk has a negative radial metallicity gradient while the thick disk has no radial metallicity 
gradient. In addition, as Fig. 9 shows, the data points in Region B are in general well described by a 
straight line of a single slope, while the mean metallicities at the inner radial annuli (< 8 kpc) show 
some deviations from this slope and tend to to exhibit a flatter slope even for height bins close to the 
midplane. Similar results have also been found recently by Hayden et al. (2014) using APOGEE red 
giants as tracers. This flattening of gradient in the inner region is possibly caused by the existence of a 
central bar that mixes stars of different populations, as suggested by Hayden et al. (2014). 

For the outer disk (i.e. Region C), as described earlier, the previous studies are not in conclusive 
with regard to whether the gradient is flatter than near the solar circle owing to the poor sampling of 
the available samples. In the current study, with a large number of RC stars sampling the outer disk, we 
are able to confirm that the metallicities in the outer Region C is essentially constant, and the gradients 
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Fig. 11 Distributions of metallicities, plotted as a function of \Z\, deduced for the individual 
annuli of Rgc from the RC sample. The red dashed lines divide the disk into two regimes: the 
lower (| Z | < 1 kpc) and upper (\Z\ > 1 kpc) disk. Red boxes gives the mean metallicities in 
the individual vertical bins. The linear fits to the data of the lower and upper disks are shown 
by green lines, while those for the whole disk height are plotted in red. The gradients yielded 
by the fits are also marked in each plot, with the color code as for the fit. 
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Fig. 12 Vertical gradients, A[Fe/H]/AZ, plotted as a function of /ice, deduced for different 
ranges of disk height \Z\ (Red: \Z\ < 3kpc; Blue: \Z\ < 1 kpc; and Mangeta: 1 < \Z\ < 
3 kpc). The lines connecting the data have been smoothed over 3 adjacent points. The grey 
triangles are measurements from Carrell et al. (2012) using SEGUE F/G/K dwarfs of 1 < 
\Z\ < 3kpc. 


derived are indeed less steep than those found for the solar circle. The different trends of radial gradients 
as a function of \Z\ between RegionB and C, as shown in Fig. 10, may suggest that these two parts of 
the disk may have quite different evolution paths, resulted from, for example, infalling gas enriched 
in the halo (e.g. Chiappini et al. 2001), the presence of non-axisymmetric structures (e.g. a central bar, 
long-lived spiral arms; e.g. Scarano & Lepine 2013), or a merge event in the outer disk (e.g. Yong 2006). 

4.2 Vertical metallicity gradients 

The vertical gradients are poorly determined in previous studies. Most of those studies find a steep 
gradient near the solar circle, for example, —O^Sdexkpc^ 1 by Bartasiute etal. (2003), —0.30dexkpc 1 
by Chen et al. (2003) and —0.29 dex kpc -1 by Marsakov & Borkova (2006). Those values are all slightly 
steeper than our result near the solar circle. The discrepancies are possibly caused by the fact that the 
sample tracers used by those earlier studies cover a wide range of Rgc and the vertical gradients are 
deduced without binning the data in Rug- Thus the values derived may be affected by the presence of 
radial gradient, the so-called the correlation effects discussed earlier. The different age distributions 
of tracers used in different studies could also be responsible for the discrepancies. For the same Rgc 
or \Z\ range, our results are in excellent agreement with the previous studies. As Fig. 12 shows, the 
vertical gradients derived at different Rgc for the upper disk from our RC sample are comparable to the 
measurements by Carrel et al. (2012) using SEGUE F/G/K dwarfs of heights 1 < \Z\ < 3 kpc. 

In the current study, we find for the first time that the vertical gradients of lower disk (0 < \Z\ < 
1 kpc) flatten with increasing Rgc quicker with Rgc than that of upper disk (1 < \Z\ < 3 kpc). The 
result should provide important constraints on the formation scenarios of thin and thick disks. 

5 CONCLUSIONS 

With improved log g measurements obtained with a KPCA method trained by accurate asteroseismic 
data from the L A MOST - Keplei• fields, we have assembled hitherto the largest clean sample of RC 
stars, selected from the LSS-GAC DR2. The sample contains over 70000 objects of intermediate ages, 
with a typical distance uncertainty of 5 -10 %. The sample stars cover a significant disk volume of pro- 
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jected Galactocentric radii 7 < Rgc < 14 kpc and heights from the Galactic midplane 0 < \Z\ < 3 kpc. 
With this RC sample, we have measured the radial and vertical metallicity gradients and studied their 
spatial variations. The main conclusions of the current study are as follows: 

1. Both the radial and vertical gradients are negative across much of the disk probed. The gradients 
also show significant spatial variations; 

2 . Near the solar circle (7 < Rgc f£ 11.5 kpc), the radial gradients flattens as one moves away from 
the Galactic plane; 

3. In the outer disk (11.5 < Rgc < 14kpc), the radial metallicity gradients do not show clear varia¬ 
tions with \Z\, and have an essentially constant value of —0.014 dexkpc ' 1 much shallower than near 
the solar circle. This suggests that the outer disk may have a very different evolution path; 

4. The vertical gradients also show significant spatial variations, flattening with increasing Rgc- 
More interestingly, for the first time, we find the vertical gradients of the lower disk (0 < \Z\ < 1 kpc) 
flattens with -Rgc quicker than those of the upper disk (1 < \Z\ < 3 kpc). 

For the moment, the chemical or chemo-dynamical models can not fully predict the spatial vari¬ 
ations of metallicity gradients found in the current work. Our results suggest that processes such as 
radial migration, non-axisymmetric perturbation, and galaxy mergers may have played a significant role 
in shaping the Galactic disk(s). The LSS-GAC survey is an ongoing project and will last until 2017. 
With more observations becoming available, the size and spatial coverage of our sample will continue 
to expand. In addition, additional information, especially the [a/Fe] ratios, will become available with 
the new version of LSP3. As the survey progresses, the LSS-GAC survey will in no doubt continue to 
provide further, stronger constraints on the chemical enrichment history of the Galactic disk(s), and shed 
light on the formation and evolution of the Galaxy. 
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Appendix A: 

The number of PCs adopted in our current KPCA model, i.e. Npc = 25, is a tradeoff of the following 
two considerations. Firstly, a sufficient number of PCs is required to construct a tight relation between 
asteroseismic logy and LAMOST blue-arm spectral features. As Fig. Al shows, when Npc is small, 
the fit residuals show a large systematics. When JVpc > 25, the fit residuals no longer show obvious 
trend as the atmospheric parameters vary and the standard deviations of the fit residuals decrease slowly 
with increasing Npc- On the other hand, a value of Npc that is too large will lead to biased estimates of 
log g under low spectral S/N’s. As Fig. A2 shows, the systematics become quite significant for spectral 
S/N’s smaller than 20 if Npc > 30. Hence, we have adopted in the current work Npc = 25 for 
KPCA model in constructing the relation between asteroseismic log g and LAMOST blue-arm spectral 
features. 


This paper was prepared with the RAA I4TpX macro vl.2. 
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Fig.A.l Distributions of values of residual, logpKPCA—loggAST, of the training sample, as 
a function of LSP3 T e g and [Fe/H] and of asteroseismic log g for different iVpc assumed. The 
last panel in each subplot shows a histogram distribution of the residuals (black line). Also 
overplotted in red is a Gaussian fit to the distribution. The mean and dispersion of the fit, as 
well as the number of stars used, are marked. 
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Fig. A.2 Differences of KPCA log g estimates and asteroseismic values for the control sam¬ 
ple, plotted as a function of spectral S/N (4650 A) for different Npc assumed. The blue dashed 
lines indicate differences A log g = ± 0.1 dex. The means and standard deviations of the dif¬ 
ferences in the individual S/N (4650 A) bins (binsize = 8 ) are oveplotted with dots and error 
bars, respectively. 























